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ABSTRACT 

A Kerr black hole can flare quasi-periodically, if it is magnetically connected to 
the accretion disk. The flares extract the rotational energy of the hole. The flaring 
occurs because the hole twists the magnetic field lines beyond the threshold of the 
screw instability. 

Subject headings: black holes - magnetic fields 

1. Introduction 

We show that a Kerr black hole can flare quasi-periodically, if it is connected to the accretion 
disk by a magnetic fleld. The averaged luminosity in the flares is ~ B'^r^c, where B is the magnetic 
field at the event horizon, r+ is the event horizon radius, c is the speed of light. 

The mechanism considered here is different from the Blandford-Znajek process - the magnetic 
field connects the hole to the accretion disk, instead of to infinity. In §2 we give the astrophysical 
formulation of the problem and an order-of-magnitude solution. In §3 we idealize the problem and 
solve it rigorously. 

2. Astrophysics of Black Hole Flares 

Fig. 1. shows a Kerr black hole, its accretion disk, and a magnetic spot in the disk. The hole 
will twist the magnetic field lines, that is a toroidal magnetic field will be generated. To calculate the 
toroidal field, we treat the hole as a rotating conducting ball, radius ~ r^, conductivity o" ~ c/r+ 
(Znajek 1978). We assume for simplicity that the ball is rotating at about the speed of light, 
corresponding to a nearly maximally rotating black hole. In the ball rest frame, the poloidal 
electric field is E ^ B, where B is the magnetic field. The current is / ~ r'^aE ~ cr^B. This 
current creates a toroidal magnetic field, Bt ~ I/{cr^) ~ B. 

The twisted magnetic field will try to unwind itself. According to the Kruskal-Shafranov 
criterion (Kadomtsev 1966), this leads to the screw instability when the toroidal field becomes so 
strong that the magnetic field line turns around itself about once. The instability condition is 
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Fig. 1. — Kerr black hole in a magnetic contact with the accretion disk. 

Bt/vj^ ^ B/L^ or L ^ r^, where L is the distance from the event horizon to the magnetic spot 
measured along the field line. 

If the stationary magnetic configuration is screw-unstable, the black hole magnetosphere fiares 
and releases part of the magnetic energy. Then the hole twists the magnetic field again, and 
the magnetosphere becomes screw-unstable. The flares repeat quasi-periodically, with the period 
^ r+/c. 

3. Idealized problem 

We straighten the field lines connecting the hole and the magnetic spot in the disk. We neglect 
the motion of the spot as compared to the hole rotation. This transforms Fig.l into Fig.2. 

Consider a Kerr black hole and two ideally conducting plates at a distance L from the hole. 
The plates are perpendicular to the black hole angular momentum. There is a uniform external 
magnetic field Bq frozen into the plates. The electromagnetic field is force-free. We calculate 
the force-free magnetosphere in §3.1, formulate the stability problem in §3.2, and show that the 
magnetosphere can be screw-unstable in §3.3. 

The results given here constitute a rigorous proof that the magnetosphere of a Kerr black hole 
is screw-unstable for certain values of L ^ ^ a, where a is the angular momentum of the hole 
per unit mass. In fact, this means that a fast, a ~ r+, Kerr hole will flare if L > Cr_)_ where C is a 
dimensionless number of order unity. 
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Fig. 2. — Kerr black hole in a magnetic field. At a distance L from the black hole, the magnetic 
field is frozen into ideally conducting plates. 



3.1. The force- free magnetosphere 

Wc derive the equation for the force- free magnetosphere following the steps similar to Blandford 
&; Znajek (1977). In our case the equations and their derivation are much simpler because the 
electrostatic potential is everywhere zero due the ideally-conducting plates. 

The Boyer-Lindquist coordinates are used, with c = 1, and G = 1. The metric is 



2Mr 2 4Mar . 2 



ds^ = (1 - -^)dt^ + 



^ sin^ edtd^ - ^dr^ - J^de^ - + ^ ^ sin^ Od^^ (1) 



S = + cos^ 9, 

A = - 2Mr + a^, 



(2) 
(3) 



(4) 



M is the mass, a is the angular momentum per unit mass. The event horizon radius is r+ 
M + VM2 - a2. 

The magnetosphere is described by the Maxwell equations 

dkiV^F"') = -47rv^/, 

and the force-free condition 

Fikj^ = 0. (5) 

Here g = det{gik) = — S^sin^^, Fik = diAk — d^Ai is the electromagnetic field tensor, A^ is the 
four-potential, f is the four-current. We consider the stationary axisymmetric case, dt = = 0. 

Prom (5), 

{fdr+i'dg)At = 0, (6) 
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{fdr+fde)A^ = f). (7) 

Then At must be a function of A^. At the plates At = 0. Assuming that aU the points in 
the magnetosphere are connected to the conducting plates by the magnetic field lines, = 
everywhere. 

Denote A = —A^, B = —drAg + deAr. Then the remaining two force- free conditions and three 
of the Maxwell equations are 

j't'drA + Bf = 0, (8) 

fdeA - Bf = 0, (9) 

de{V^g"'9''B) = 47^^/^ (10) 

driV^g'-'g^'B) = -4v^^J^ (11) 

driV^g^^g'-'drA) + deiV^g'^'^'g^'deA) = -AnV^jt (12) 

The combination ^/—gg^^g^^B = 21, where / is the poloidal current. From (10) and (11), 

{fdr+}'de)I = 0. (13) 

From (7) and (13), / = I (A) - the current flows along the field lines. From (10), —gf = I'dgA, 
where /' = dl/dA. Now from (9), 

TTggrrgm 

With j'^ given by (14), j*" by (10), and by (11), equations (7)-(ll) are satisfled. Equation 
(12) takes the final form 

drig^odrA) + smede{{gooh^^O)deA) + AA-'^II' = 0, (15) 

where 500 = 1 ~ 2MrS^^. The boundary conditions at the plates are 

AU=±L = Bop^/2, (16) 

where z = rcosO, p = rsinO. The boundary condition at the event horizon follows from the 
requirement that F'^ be finite at r+. This gives 

2I=^dgA\r=r^. (17) 

We now solve (15) assuming a <C r+. For a = 0, (17) gives / = 0, and then (15) and (16) give 
A = Bqt^ s\v? 6/2. We use this A in (17) and find the poloidal current. 



, clA 2A , ^ 

I{A) = —Jl-—^. 18 



We can now use this I{A) in (15) and iterate if greater precision is required. But we do not need 
a more accurate solution for our current purpose. Equation (18) is sufficiently accurate to analyze 
the stability of the magnetosphere in the limit of r+ ^ a. 



^The remaining Maxwell equation gives the charge density. We do not need it. 
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3.2. Stability in Force-Free Electrodynamics 

We have found the stationary magnetic configuration. Now we study its stability in the frame- 
work of force- free electrodynamics (FFE) defined by (4), (5). For L ^ r+, we can use the galilean 
metric. 

FFE is classical electrodynamics supplemented by the force-free condition. Thus 

dtB = -V X E, (19) 

9tE = VxB-j, (20) 

/9E + jxB = 0. (21) 

V • B = is the initial condition. The speed of light isc=l;p = V- E and j are the charge and 
current densities multiplied by 47r. The electric field is everywhere perpendicular to the magnetic 
field, E • B = 0. The electric field component parallel to the magnetic field should vanish because 
the charges are freely available in FFE. It is also assumed that the electric field is everywhere 
weaker than the magnetic field, E'^ < B^. Then equation (21) means that it is always possible to 
find a local reference frame where the field is a pure magnetic field, and the current is fiowing along 
this field. FFE is Lorentz invariant. 



Equation (21) can be written in the form of the non-linear Ohm's law 
._ (B-VxB-E-Vx E)B + (V • E)E x B 

J - ^2 • 



(22) 



Equations (19), (20), (22) form an evolutionary system (initial conditions V • B = and E • B = 
are implied). It therefore makes sense to study stability of equilibrium electromagnetic fields in 
FFE. One can also study linear waves and their nonlinear interactions in the framework of FFE 
(Thompson & Blaes 1998). 

It is convenient to introduce a formulation of FFE similar to magnetohydrodynamics (MHD); 
then we can use the familiar techniques of MHD to test stability of magnetic configurations. To 
this end, define a field v = E x B/S^, which is similar to velocity in MHD. Then E = — v x B and 
equation (19) becomes the "frozen-in" law 

OjB = V X (v X B). (23) 

From V = E X B/S^, and from equations (19)-(21), one obtains the momentum equation 

dtiB"^^) = VxBxB-hVxExE+(V- E)E, (24) 

where E = — v x B. Equations (23), (24) are the usual MHD equations except that the density is 
equal and there are order corrections in the momentum equation. 



^If p 7^ 0, E ■ B = follows from (21), if p = 0, this condition is an independent basic equation of FFE. 
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A force-free magnetic field (V x B x B = 0) with a zero electric field is a stationary FFE 
solution. We will study the stability of the force-free magnetic configurations using the MHD 
formulation of FFE. Linear perturbations, denoted b and v, oc exp(— iwi), satisfy 

- zwb = V X (v X B), (25) 

- iujB\ = VxBxb-FVxbxB. (26) 

Define the displacement ^ by v = — io;^. From (25), b = V x (^ x B). Now equation (26) can be 
written as 

- u^B'^$, = VxBxVx(|xB)-FVxVx(|xB)xB = -K^. (27) 

Since the operator K is self-adjoint (Kadomtsev 1966 and references therein), the frequency is 
given by the variational principle 

The equilibrium field B is unstable if the potential energy 

W = J d^r^K^, (29) 

is negative for some displacements ^. 

Consider the case of cylindrical symmetry, coordinates (r,9,z). The equilibrium field B = 
(0, U{r), B{r)) should satisfy BB' + Ur-^{rU)' = 0, where the prime denotes the r-derivative. For 
an eigenmode oc exp(im0 -|- ikz), the potential energy reduces to (Kadomtsev 1966) 

j drue+ge), (30) 

where ^ is the radial component of the displacement. The other two components of the displacement 
vector, and (^z^ were chosen to minimize the energy for a given ^. The functions of radius / and 
g are given by 

j_^ {krB + mUf ^ 

^2^2 _j_ j^2 ' 

We now show that magnetic configurations with a non-zero toroidal field U are screw unstable. 
Screw means, e. g., that m = 1, but not m = —1, is unstable for a given sign of k. Kadomtsev 
(1966) gives a clear discussion of the screw instability in plasmas, and shows that the screw mode 
is the most dangerous mode (if the plasma is unstable, it is screw unstable) . 

Assume that C/(0) = 0, C/'(0) > ^, U — for r — oo, and U is positive in between. Assume 
that B is everywhere positive. Let k be positive and small. Take m = — 1. Let tq be the first 
zero of krB — U. Take a trial function ^ = 1 for r < ro, and zero at r > tq. We can chose this 
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generalized function ^ in such a way that the first term in the energy integral (31) vanishes. The 
second term is an integral of g from to tq. For small k, g = k^r{krB — U){3krB + U) is negative. 

Thus an eigenmode a exp(— + ikz) is screw-unstable if krB ^ U. If the length of the 
magnetosphere in the z-direction is L, the minimal value of A; is ~ L~^. The instability condition 
(the Kruskal-Shafranov criterion) is 

- > C^, (33) 
r L 

where C is a dimensionless number of order unity. The value of C depends on the boundary 
conditions. 



3.3. Stability of the Black Hole Magnetosphere 

For the black hole magnetosphere, in the limit r+ S> a, at a distance from the hole S> r+, the 
poloidal field B is approximately constant and equal to Bq. The poloidal current is given by (18), 
and the toroidal field is 




(34) 

in the cylindrical coordinates of §3.2. According to the Kruskal-Shafranov criterion, the magneto- 
sphere will be screw unstable if L > Cr^/a, where C is a dimensionless number of order unity. 



4. Discussion 

The flaring is a new mechanism for the extraction of the rotational energy of Kerr holes. Many 
previous studies considered the Blandford-Znajek mechanism in which the magnetic field connects 
the event horizon to infinity. By contracts, the fiaring occurs when the magnetic field ends in the 
accretion disk. Both mechanisms give the same averaged luminosity, as expected from dimensional 
analysis. 

The essential feature of our mechanism is that the rotational frequency of the hole determines 
the time structure of the radiated energy. The hole frequency may appear observationally as a 
quasi-periodicity of the X-ray flux. It would be interesting to search for this periodicity in the 
high-time-resolution data. 
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